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We consider an ensemble of three-level particles in lambda-configuration interacting with two 
bosonic modes. The Hamiltonian has the form of a generalized Dicke-model. We show that in the 
thermodynamic limit this model supports a superradiant quantum phase transition. Remarkably, 
this can be both a first and a second order phase transition. A connection of the phase diagram to 
the symmetries of the Hamiltonian is also given. In addition, we show that this model can describe 
atoms interacting with an electromagnetic field in which the microscopic Hamiltonian includes a 
diamagnetic contribution. Even though the parameters of the atomic system respect the Thomas- 
Reiche-Kuhn sum rule, the system still shows a superradiant phase transition. 
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I. INTRODUCTION 

Superradiant phase transitions in systems of atoms 
interacting with an electromagnetic field were first dis- 
cussed theoretically by Hepp and Lieb [TJd], and by Wang 
and Hioe (3] . The superradiant phase is characterized by 
a macroscopic and coherent excitation of both the atoms, 
and the electromagnetic field modes. To date, superradi- 
ant phase transitions have been observed experimentally 
in artificial realizations of the seminal Dicke-model |3] in 
cold atoms |5J, but not for real or for artificial atoms. The 
question which was debated forty years ago and which is 
under debate today again, is whether or not superradi- 
ant phase transitions are in principle possible in these 
systems [SHE]. 

The system of an ensemble of atoms interacting with 
an electromagnetic field is often described by the Dicke- 
model [4 . In this model, each individual atom is de- 
scribed by a two-level system and the electromagnetic 
field by a single mode of a resonator. For the general- 
ized Dicke-model with diamagnetic terms included (i.e. a 
Hopfield-like model [3] ) , there exists a no-go theorem 
|8] which precludes the transition to a superradiant phase. 
This no-go theorem is a consequence of sum rules and of 
the appearance of a diamagnetic term in the microscopic 
Hamiltonian which is quadratic in the transverse vector 
potential A. Thus no matter how low the temperature 
or how strong the coupling strength is, the superradiant 
phase transition does not exist in such models. 

The no-go theorem has also been extended to multi- 
level Hopfield- like-models j8|. However, it has been 
commented [TU] that these kind of no-go theorems for 
multi-level systems do not apply to first order superra- 
diant phase transitions and therefore new perspectives 
open up. 

In this article we extend the abstract model of our pre- 
vious paper |llj by including a diamagnetic term and ac- 
count for all boson mediated transitions from the ground 
states to the excited state. We analyze and discuss the 
phases and phase transitions of this extended model, and 
show that the system shows both a first and a second or- 




Figure 1: (Color online) Lambda-configuration of energy lev- 
els of a single particle: The two bosonic modes with energies 
uj\ and L02 couple the excited state with energy E$ to either 
of the two ground-states with energies E\ and E2 of the par- 
ticles, respectively. The corresponding coupling strengths are 
given by g\ and g2- The geometry of the setup where the 
particles are atoms and the bosonic modes are modes of an 
electromagnetic field (cf. Sec. |A| is shown on the top. In our 
numerical analysis we assume that the polarization vectors e„ 
are parallel to the dipole matrix elements d3„, respectively. 
In addition, the polarization vectors must not be orthogonal. 



der phase transition. We also discuss the symmetries of 
this model and show how they are related to the quan- 
tum phases of this system. Finally, we map this model to 
a system consisting of atoms interacting with an electro- 
magnetic field and show that a no-go theorem in gen- 
eral does not hold for this model and a superradiant 
phase transition is in fact possible. However, due to 
the Thomas-Rcichc-Kuhn sum rule, only the first or- 
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der phase transition survives. This is consistent with the 
results of Ref. |12) . There, three-level atoms were con- 
sidered as well, and only a first order superradiant phase 
transition was found. 



II. THE MODEL 

We consider a system consisting of N identical parti- 
cles interacting with two modes of a bosonic field. The 
particles are described by three-level systems in lambda- 
configuration, i.e. two in general non-degenerate ground- 
states and one excited state with energies E\, E 2 , and E 3 , 
and detunings 5 = E 2 — E\ > and A = E 3 — E\ > 0. 
The two bosonic modes have energies oj\ and uj 2 , respec- 
tively. In Fig. [I] the model is summarized graphically. 
The Hamiltonian is given by (h = 1) 
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Here, J n m are collective atomic operators (see appendix |a| 
or Ref. and a\, a n are the creation and annihilation 
operators of the nth bosonic mode. The first two terms 
in Eq. ([IJ correspond to the energy of the free system of 
particles and a bosonic field, respectively. The third and 
the fourth term proportional to the coupling constants g\ 
and g 2 model the interaction between the particles and 
the bosonic field. The dimensionless parameters xi (X2) 
generate transitions from the first (second) ground state 
to the excited state induced by the second (first) bosonic 
mode. Eventually, the last terms scaling with K\, k 2 , and 
k 3 represent a self interaction of the bosonic field. 

Special cases of the Hamiltonian, Eq. ([!]), include (i) 
the model of our previous publication with Xn = 
K n = 0, and (ii) a model describing atoms interacting 
with two modes of an electromagnetic field with 
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The latter mapping and the detailed definition of the 
parameters is given in appendix [5] 

By virtue of the Thomas -Reiche-Kuhn (TRK) sum 
rule, the couplings g n are bounded. In our model, we 
have two important sum rules (the derivation can be 
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Thus for given parameters, the light-matter coupling 
strengths g\, g 2 for atomic systems can not exceed 
<?i,trk, 92, trk, respectively. This bound will be crucial 
in the next section where we discuss phase transitions of 
this model. 



III. PHASE TRANSITIONS 

In order to obtain the phase diagram of this system, 
we apply the techniques presented in our previous publi- 
cation There, we have used the Holstein-Primakoff- 
transformation [13, 14J. This is a non-linear transforma- 
tion, which maps the collective particle operators j„ m 
(which fulfill the algebra of generators of the unitary 
group) onto bosonic creation and annihilation operators. 
The corresponding non-linear Hamiltonian can be lin- 
earized by introducing mean-fields both for the bosonic 
modes, and for the creation and annihilation operators 
of the particles. Then, the Hamiltonian obtains the form 

H = NhW +VWtt 1 '> +U 2) +0{N- 1 / 2 ), (6) 

where the order of N is explicitly given. Here, corre- 
sponds to the ground-state energy, and h^> to the fluctu- 
ations around the ground-state. In the thermodynamic 
limit, N 3> 1, analysis of the ground-state energy 
yields the relevant information for the phases, phase tran- 
sitions, and phase diagram of the system. Explicitly, the 
ground-state energy is given by 

E := h<M = ^+A*|+L+^^+(wa+^)^ 
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The real mean-fields ip s (s — 1,2), $f n (n = 2,3), 
and ip\ = \f\ — ^\ — ^3 correspond to the two bosonic 
modes and the bosons introduced by the Holstein- 
Primakoff-transformation, respectively (cf. Ref. |11]). 
Finite values of the mean-fields give macroscopic pop- 
ulations of the bosonic modes or the three particle levels, 
i.e. if $3 ((fx) is finite, then the third (first) energy level 
of the particles (bosonic mode) is macroscopically occu- 
pied. Particle number conservation implies ^| + ^3 < !• 
We eliminate the two mean-fields ip\ and (p 2 from E 
by minimizing Eq with respect to these two mean-fields. 
Then, the ground-state energy is a function of the mean- 
fields ^2 and ^3 only: Eq = £0(^2, ^3)- One can show 
that the normal state with ty 2 = fy 3 = is always a 
critical point of Eq. However, it can still be a maximum, 
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or. if it is a minimum, it can be a local minimum only. 
To check whether it is a minimum or a maximum, we 
analyze the Hessian of E (^ 2 = 0, ^3 = 0). This gives 
the inequality 
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As long as this inequality is fulfilled, the normal state 
minimizes the ground-state energy. 

Analyzing the Hessian of £0(^2 = 0, ^3 = 0) we can 
check whether or not the normal state minimizes the 
ground-state energy. However, we do not see if it is a 
global minimum of Eq- To see this, we need to check 
all critical points of Eq. We do this by numerical mini- 
mization of the ground-state energy E . Fig. [2] shows the 
order parameters ^2 and ^3 obtained by this numerical 
procedure. Regions with ^ = and \&3 = correspond 
to the normal phase, whereas the superradiant phase is 
characterized by regions with f 2 > or $ 3 > 0. As is 
clear from Fig. [2] the system supports a normal and a 
superradiant phase. 

If two thermodynamic phases are separated by a sec- 
ond order phase transition the order parameter (which is 
given by the mean-fields) characterizing this phase tran- 
sition is continuous across the phase boundary. However, 
numerical analysis indicates that it is not a second order 
phase transition along the entire phase boundary, and at 
a certain point the order parameters, e.g. ^2, become 
discontinuous. This behavior is shown in the lower panel 
of Fig. [2] Then, the phase transition is of first order. 
Consequently, the phase boundary separating the nor- 
mal from the superradiant phase consists of two segments 
which characterize a first (solid green line in Fig. |2| and 
a second (dashed green line in Fig. [2]) order phase tran- 
sition, respectively. These two segments meet in a single 
point (red ring in Fig. [2|. The location of this point can- 
not be obtained by an analysis of stability of the normal 
phase. 

The continuous phase transition signals the breakdown 
of stability of the normal phase. Therefore, the phase 
boundary for the second order phase transition is given 
by the value of g\ where the inequality, Eq. (|8| , becomes 
an equation, i.e. for g\ — <? lc . Since g\_ c does not depend 
on <?2, the phase boundary is a straight line. With the nu- 
merical values used to generate Fig. [2j gi tC — 1.03 (^trk 
(dashed red line in Fig. |2|, where <?i,trk is given by 
Eq. (5). This is in perfect agreement with our numer- 
ical findings (cf. Fig. |2j. 

Concerning the atomic system, the couplings g n are not 
arbitrary and can not exceed a bound given by the TRK 
sum rule of Eqs. These bounds are displayed by the 
solid yellow lines in Fig. [2j As can be seen in Fig. [2] this 
atomic system can in fact undergo a superradiant phase 
transition depending on the choice of the parameters. 
The numerical calculations show that this is a first order 
phase transition. Indeed, combining the TRK sum rule 



for gi, Eq. ^yj, with the stability criterion, Eq. (|8j> , for 
the normal state, we see that in our model the normal 
state is stable for all gi < ^i^trk- Hence, concerning our 
model, the superradiant phase transition can not be of 
second order. This is also in agreement with the results 
found for the single-mode model of Ref. [12_. 



IV. SYMMETRIES 

If one neglects the terms proportional to n n (n — 
1,2,3) and Xn ( n = E2), then the Hamiltonian, 
Eq. ([T]), commutes with the symmetry-operators = 
exp[— i7r(— i n ™ +a^a„)], (n — 1,2). As was shown in 
Ref. [TT] , this gives rise to two superradiant states where 
a single bosonic mode is macroscopically occupied only 
(blue and red superradiant states), and either of the 
superradiant phases correspond to one broken symme- 
try. In contrast, including those terms, the symmetry- 
operators do not commute with the Hamiltonian, Eq. (JlJ, 
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Hence, the parities corresponding to the two branches of 
the lambda-system are no conserved quantities anymore, 
since the term proportional to (di + a\)(a,2 + a\) 'mixes' 
both bosonic modes. However, the operators 

fl„ = exp[-wr(-7 n " + a\ ai + a\a 2 )\ (n = l,2) (10) 

do commute with the Hamiltonian, Eq. ([T|), and therefore 
this parity including both bosonic modes is conserved. 

In the normal phase both parities Hi, n2 are con- 
served, whereas both parities are simultaneously broken 
in the superradiant phase. This loss of parity symmetry 
has a large impact on the phase diagram: In our nu- 
merical analysis, we obtain only one superradiant phase 
with both bosonic modes being macroscopically excited 
simultaneously (cf. Fig. [2]). In addition, analytical cal- 
culations show that blue and red superradiant states are 
not stable. This is in contrast to the model without the 
terms proportional to k„ and Xn, where these superradi- 
ant states are the only superradiant phases and a state 
with both branches of the lambda-system being simulta- 
neously superradiant is not stable [11] . 



V. CONCLUSION 

In this work, we have presented a generalized Dicke- 
model of particles with three energy-levels in lambda- 
configuration coupled to two bosonic modes where the 
microscopic Hamiltonian contains a diamagnetic term. 
We showed that this system exhibits a superradiant 
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quantum phase transition in the thermodynamic limit. 
This phase transition can be both of first and of second 
order and we analytically derived the critical coupling 
strength for the second order phase transition. Quan- 
titatively, the whole phase diagram was obtained using 
numerical methods. Compared to the model without the 
diamagnetic term we have studied in a previous publica- 
tion the phase diagram has one superradiant phase 
only. The loss of the second superradiant phase is directly 
connected to the diamagnetic term since it changes the 
parity symmetry of the Hamiltonian. 

In addition, we mapped this abstract model to an 
atomic system interacting with two photonic modes of 
a resonator and showed by numerical calculations that 
the superradiant phase transition persists. We emphasize 
that the microscopic Hamiltonian includes diamagnetic 
contributions and that this model respects the Thomas- 
Reiche-Kuhn sum rule which gives bounds for the cou- 
pling strengths. This is in stark contrast to Hopfield- 
like-models, where the combination of the diamagnetic 
contribution and the sum rule suppresses the superra- 
diant phase transition. Thus, compared to the Dickc- 
model with diamagnetic terms included and its general- 
izations [5], no no-go-theorem exists in our model. How- 
ever, we showed that the superradiant phase transition 
is of first order. 

In experiments, our model would be realized by atoms, 
if they can be reduced to three-level systems in lambda- 
configuration. In addition, the dipole matrix element d 31 
(d32) must not be orthogonal to the polarization vector 
£2 (£1) of the two modes of the resonator. 

In the paper of Baksic et al. [12_ a similar model 
is discussed. They consider atoms in a general three- 
level-configuration, which are coupled to one mode of a 
resonator. Including a diamagnetic contribution in the 
Hamiltonian and respecting the Thomas-Reiche-Kuhn 
sum rule, they find a superradiant quantum phase tran- 
sition. Thus, no no-go-theorem exists in their model 
cither. In particular for the lambda-configuration, this 
phase transition is always of first order. This agrees with 
our results. 

There is one point we want to remark: Theoretically, 
the Dicke-model (i.e. a Hopfield-like-model without the 
diamagnetic term) supports a second order phase tran- 
sition [TH51 IT51 115]. The no-go theorem [6-8 applies to 
these continuous phase transitions only (see appendix [Cj. 
However, first order phase transitions could in principle 
still provide a superradiant phase. 
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Appendix A: Derivation of the model 

In the Coulomb-gauge, the Hamiltonian of N identical 
atoms interacting with an electromagnetic field inside a 
resonator is given by (cf. e.g. [T7|) (H = 1) 



H = 



N / N e 

i=i \ j=i 



Pij ~ 9 A ( r y) 



+ V c ({r l i,---,r l N c }) + 



^cj s a|a s . (Al) 



Here, m, q, p^- and are the mass, the charge, the 
kinetic momentum and the position of the jth electron 
of ith atom respectively, N e is the number of electrons 
per atom, A is the transverse vector potential, V c is the 
Coulomb-potential of all electrons with respect to their 
respective nuclei, uj s is the dispersion of the sth mode 
of the transverse electric field and a\, a s create, annihi- 
late one photon in the sth mode. In the dipole approx- 
imation, i.e. when the wavelength of the electric field 
is large compared to the dimension of the system, the 
transverse vector potential is given by (cf. e.g. |17| ) 
A = ^ s A s e s (oj + a s ) . The quantum number s contains 
both polarization of the transverse electric field and the 
momentum k of the photon. Real polarization vectors e s 
are considered only. The parameter A s is defined by 
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with the So the vacuum permittivity and V the volume 
of the resonator. 



The Hamiltonian, Eq. ( Al I, can be written as 
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Here, the eigensystem {E n , \n) 3 } of the free system - 
that is the kinetic energy of the N e electrons of a sin- 
gle atom, plus the Coulomb-energy of the N e electrons 
with its respective nuclei — has been introduced. The 
energies E n are the same for every atom. The next 
step is to express the kinetic momentum operator p^ 

in this basis. The identity p,y = im\Jii\ fy] gives 

Pij = imJ2„,i(E„ - Ei) {n\hj\l) \n) 3 (l\. 
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Next, we introduce the coupling constants 

g n i jS = -iVN(E n -Ei)A a e s -d nh (A7) 



is given by 



with d„, = (n\r l:j \l) , 



(A8) 
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where the dipole matrix elements d„; are identical for all 



atoms. Then, the Hamiltonian, Eq. (A3 1, assumes the 
form 



S 



This holds for all quantum numbers I . This kind of sum 
rule for the dipole matrix clement d n i is called Thomas- 
Reiche-Kuhn-sum rule |18j . 

Since in our model the infinite dimensional Hilbert- 
space of a single atom is restricted to three energy levels 
only, and di2 = is assumed, the left side of Eq. (B4| is 
bounded, i.e. 
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Eventually, using the definitions for g n (Eq. (All l) and 



with collective operators /„' = 2,=i \ n Y(A-< an d the pa- 
rameter k is given by 



q 2 NN e UJ n 



2m 
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Finally, we restrict the Hamiltonian, Eq. (A9l, to a 



three-level lambda-system with two modes of a resonator, 
i.e. n — 1,2,3 and s = 1,2. We assume that the dipole 
matrix element d 12 = 0. This may be due to symmetry. 
Eventually, introducing real coupling constants 



for k (Eq. (AlOl), and choosing I = 1,2, we obtain the 
two inequalities of the Eqs. 



Appendix C: No-go theorem for second order 
superradiant phase transitions 

We consider the most general model of an ensemble of 
N multi-level atoms interacting with a single mode of a 
resonator. The Hamiltonian has the form 



g n = VN(E 3 -E n )\d 3n -e n \A n , (n=l,2) (All) H=J2 E Jr 



" + uoa) a + — (a T + af 



it J. ^2 



the Hamiltonian of Eq. with Xn = ^ (n/! = 

1, 2), a n i = |d 3n • £;|/|d 3 „|, ki = k 2 — k, k 3 — Kyfa, and 
a = Ei ■ £2 is obtained. 



Appendix B: Derivation of the 
Thomas Reiche Kuhn sum rules 

For every Hamiltonian H with spectrum E n and eigen- 
basis \n) the identity 

J2(En - El ) (l\d\n) (n\6\l) = \ (l\ [6, [H, 6}] \l) (Bl) 



is fulfilled. Here, the sum runs over all quantum numbers 
n. This equation is valid for any Operator O and every 
quantum number I of the Hamiltonian H. 

A special sum rule for atomic systems with a Hamilto- 
nian of the form 

H = Jl^a+H{^---^}) (B2) 

3=1 
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(/ n w +/ m «)(at+a). (CI) 



Each single atom has v energy levels characterized by the 
non-degenerate energies E\ < E2 < ■ ■ ■ < E v . Popula- 
tions of and transitions among the atomic energy levels 
are described by the collective operators I n n and / n m , 
respectively. The energy of the mode of the resonator 
is given by w > 0; the bosonic operator (a) create 
(annihilate) a corresponding photon. The coupling of 
the transverse vector potential with the atoms is approx- 
imated by the dipole coupling, with coupling strength 
9n,m — 3m. n, 9n,n = 0. The diamagnetic contribu- 
tion is parameterized via the parameter k. In deriv- 
ing this Hamiltonian, we have neglected atom-atom and 
Coulomb- interaction among different atoms. 

The TRK sum rules (cf. Eq. (B4|) of this model can 
be written in the form 



a* 



, E n 

n=l 
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and operators 



,V. 



6 = e-d, d = qJ2^i, N=l (B3) 



with E„, m = E n — E m and m = 1, . 

2 

terms of oscillator strengths f n , m = g "' r " K 2 
sum rule obtains the memorable form ^ n /„ 



(C2) 



Note that in 
^, the TRK 
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Applying a Holstein-Primakoff-transformation to the 
collective operators I n m and introducing mean fields both 
for the atomic [f n , (n = 2, . . . , v)\ and for the photonic 
(ip) degrees of freedom, we obtain the ground-state en- 
ergy per atom in the thermodynamic limit 



E, 



GS 



Eh. 



ra=2 
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E 
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Since we are only interested whether or not the normal 
state with ty„ = p = minimizes the ground-state en- 



ergy Eq. ( C3 1 , we expand the ground-state energy to 



second in order in ty n and tp. This yields 

v / ..2 \ 

E GS = E x + E n,i^l + (w + 4- 

~Vi¥ n ¥> + 0(*£). (C4) 
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Because the ground-state energy has no terms linear in 
or p. the normal state represents a critical point. In 
order to specify the type of the critical point, we calculate 
the Hessian of the ground-state energy, Eq. (C4|. We 
obtain 



^gs 



dp 2 
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d 2 E t 
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d 2 Er 



dp>df 7 



= 2^ + 4-J+0(^), 
= 2E n , x 5 n , m + 0(*„), and 
= 4g n>1 + <D(* 2 n ) 
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(C7) 



with 1 < n, m < v. Hence, the Hessian for the normal 
state is given by 



H(E GS ^ n = 0,<p = 0) 
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Now we will proof that H is in fact positive definite which 
means that the critical point is a local minimum. There- 
fore we compute the principal minors of H. The kth 
principal minor of a n x n matrix M is the determinant 
of the matrix where the first n — k rows and columns of 
M are deleted. If all principal minors of a symmetric 



matrix are positive, then the matrix is positive definite. 
In the following we will show that all principal minors of 
H are positive. 

One sees readily that the first v — 1 principal minors 
are positive since E n i are positive by definition. Hence, 
we need to compute the z^th principal minor which is 
the determinant of H itself. This is done by reducing H 
to a triangular matrix using elementary row operations. 
Eventually, this gives 
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Det[H(S GS ,^„ = 0,¥> = 0)] = 
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where X = ^ + 4^-4^= 



gn.l 

2 E„,i 



Applying the TRK 



sum rule, Eq. (C2|, with m = 1, we obtain 



X > LJ 



(CIO) 



which is always positive. Since the determinant of a tri- 
angular matrix is the product of its diagonal elements, 
the determinant of the Hessian for the normal state is 
positive. Hence we have shown that all principal minors 
are positive. Consequently, the Hessian is positive def- 
inite. This means that the normal state minimizes the 
ground-state energy irrespective of the parameters of the 
Hamiltonian and no additional states, i.e. superradiant 
states, can evolve from the normal state continuously. 
Thus, no continuous, i.e. second order, phase transitions 
are possible. However, this argument does not apply to 
first order phase transitions, and we can not say whether 
or not first order phase transitions can occur. In addition, 
a second order superradiant phase transition originating 
from a superradiant phase which results from a first order 
superradiant phase transition is also not covered by our 
analysis. We also note that a superradiant phase where 
a single one-particle energy level is macroscopically oc- 
cupied only can be excluded by a similar argument as 
presented above. 

In summary, we have considered the most general 
model of multi-level atoms interacting with a single mode 
of a resonator which obey the Thomas-Reiche-Kuhn sum 
rule. We have shown that the normal state - that is 
the state where all atoms occupy their respective ground 
state and in the resonator no photon is excited - does 
always minimize the ground state energy. Hence the nor- 
mal phase is stable irrespective of the parameters of the 
system and no second order order phase transitions to 
superradiant phases are possible. 
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Figure 2: (Color online) Top: Mean-fields ^2 and $3 for pa- 
rameters 5 = 0.1, A = 1, cji = 0.5, 0J2 = 0.6, a = V0.5, 
ki = K2 = 1, H3 = and an = an = 1- The green 

line illustrates the phase boundary and is separated into two 
segments: The solid (dashed) green line corresponds to a 
first (second) order superradiant phase transition. The point 
where these two segments meet, is marked with the red ring. 
The numerical value of gi for the boundary of the second or- 
der phase transition is given by l.OSg^TRK which coincides 
with the value obtained by an analytical analysis. The blue 
(dark) region below the green line in both diagrams represents 
the normal phase (i/>i = 1, ^2 = ^3 = f>i = <P2 = 0). The 
complementary region corresponds to the superradiant phase 
where all mean-fields are finite. The accessible parameter- 
region for atoms is indicated by the solid yellow box. The 
small inset clarifies that the superradiant phase is within this 
region. Bottom: Mean-field V&2 as a function of gi for dif- 
ferent values of (72, indicating the change of the order of the 
superradiant phase transition. The value of gi for each single 
line increases in direction of the arrow from to 1.5 with an 
increment of 0.15. In the range < gi < 0.45, ^2 is continu- 
ous as a function of g\ (solid blue lines), whereas for gi > 0.6, 
$2 is discontinuous (dashed green lines). The dashed red line 
marks the critical value for g\ of Eq. Q above which the nor- 
mal phase is unstable and correspond to the red ring in the 
top panel. In both panels, the couplings g n are scaled with 
Sn,TRK- These are the largest possible couplings for the anal- 
ogous atomic model given by the TRK sum rules, Eq. jBJ. 
The diagrams for ipi and tpi look qualitatively the same as 
the diagram for $3. In addition, ip\ is obtained by using the 
relation tpi = ^/l - \P| ~ ■ 



